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1. Introduction 

With the advent of chiral perturbation theory (ChPT), the low-energy effective field theory 
(EFT) of QCD, high accuracy calculations for hadronic reactions with a controlled error estimation 
have become possible [1,2]. The framework has been successfully applied to study, in particular, 
Tin [3] and tcN [4] scattering observables as well as nuclear forces [5]. In this contribution we 
discuss an application of ChPT to the reactions involving pion production on two-nucleon systems. 
This type of reactions allows one to test predictions of ChPT in the process with the large momen- 
tum transfer typical for the production process. As was first advocated in Refs. [6,7], the initial 
nucleon momentum in the threshold kinematics sets the new "small" scale in the problem, namely, 
p ~ sfm^Mx ~ 360MeV (% ~ p/m^ ~ ^jM n jm^), where M n {m^) is the pion (nucleon) mass. 
The proper way to include this scale in the power counting was presented in Ref. [8] and imple- 
mented in Ref. [9], see Ref. [10] for a review article. As a consequence, the hierarchy of diagrams 
changes in the modified power counting scheme of Ref. [8], and some loops start to contribute al- 
ready at NLO for s-wave pion production. On the other hand, /j-wave pion production is governed 
by the tree-level diagrams up to NNLO. In what follows we discuss the status of the theory for pion 
production in the isospin conserving and isospin violating case. It was first argued in Ref. [11] that 
charge symmetry breaking (CSB) effects recently observed experimentally in pn — > dn () [12] pro- 
vide an access to the neutron-proton mass difference which is the fundamental quantity of QCD. 
To take this challenge, however, one needs to have the isospin conserving case fully under control. 
In sec.2 we briefly discuss the theoretical status for s-wave pion production. Sec. 3 is devoted to 
a more detailed discussion of the recent results for p-wave pion production. In sec. 4 we briefly 
highlight the recent developments in the study of charge symmetry breaking effects in pn — > die . 
We finalize with the summary of the latest results. 

2. s-wave pion production and the concept of reducibility 
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Figure 1: Complete set of nucleonic diagrams up to NLO. Note that sum of all loops at NLO vanishes. 

A method how to calculate processes on few nucleon systems with external probes was pro- 
posed by Weinberg [13]: 

1 . the perturbative transition (production) operators have to be calculated systematically using 
ChPT. They should consist of irreducible graphs only. 

2. the transition operators have to be convoluted with the non-perturbative AW wave functions. 



2 



Pion reactions with few-nucleon systems 



Vadim Baru 



Figure 2: Comparison of our results to experimen- 
tal data for AW — > d%. The data sets are from 
Refs. [17-21]. Note that the green diamond and the 
blue triangle correspond to the most resent measure- 
ments from pionic deuterium atom [20,21] at rj = 0. 
The dashed line corresponds to the model of Koltun 
and Reitan [22]. The solid red curve represents our re- 
sults, as given in Ref. [14], the filled red box demon- 
strates the theoretical uncertainty of the NLO calcula- 
tion. 



Therefore it is necessary to disentangle those diagrams that are part of the wave function from those 
that are part of the transition operator. In complete analogy to AW scattering, the former are called 
reducible and the latter irreducible. The distinction stems from whether the diagram shows a two- 
nucleon cut or not. Thus, in accordance to this rule, the one-loop diagrams shown in Fig. |l](b)-(d) 
are irreducible, whereas diagrams (a) seem to be reducible. This logic was used in Ref. [9] to 
single out the irreducible loops contributing at NLO. The findings of Ref. [9] were: 

• For the channel pp — » ppn Q the sum of diagrams (b)-(d) of Fig. [j] vanished due to a cancel- 
lation between individual diagrams 

• For the channel pp — > dn + the same sum gave a finite answer 1 : 

A ^ = 4f ( - 2+3+0)lfl= i| <21) 

where f % denotes the pion decay constant and gA is the axial- vector coupling of the nucleon. 

The latter amplitude grows linearly with increasing final AW-relative momentum \q\, which leads 
to a large sensitivity to the final AW wave function, once the convolution of those with the transition 
operators is evaluated. However, the problem is that such a sensitivity is not allowed in a consistent 
field theory as was stated in Ref. [15]. The solution of this problem was presented in Ref. [14] 
(see also Ref. [16]). It was pointed out in Ref. [14] that the diagrams that look formally reducible 
can acquire irreducible contributions in the presence of the energy-dependent vertices (or time- 
dependent Lagrangian densities). Specifically, the energy dependent part of the leading Weinberg- 
Tomozawa (WT) 7tN — ► 7iN vertex cancels one of the intermediate nucleon propagators (see the 
one with the red square in Fig. |I|), resulting in an additional irreducible contribution at NLO. It 
turned out that this additional irreducible contribution compensates the linear growth of diagrams 
(b)-(d) thus solving the problem. Thus, up to NLO, only the diagrams appearing at LO (see Fig. 
|l]), contribute to pp — > dn + , with the rule that the kN — * 7iN vertex is put on-shell. The latter 
stems from the observation that in addition to the leading WT-term (~ 3M^/2) the nucleon recoil 
correction to the WT 7tN vertex (~ M K /2) also contributes at LO. As a result the dominant 7tN- 
rescattering amplitude is enhanced by a factor of 4/3 as compared to the traditionally used value 
(~ 3M K /2), which leads to a good description of the experimental data for pp — > dn + (see Fig. |2|). 

'The connection of the amplitude A to the observables is given, e.g., in Ref [14] 
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Figure 3: Diagrams that contribute to the p-wave amplitudes of AW — > NNk up to NNLO. 
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Figure 4: Results for the analyzing power at tj=0.14 (left panel) and the analyzing power at 90 degrees 
(right panel) for the reaction pp — > dn + for different values of the LEC d (in units 1 / (f^M^)) of the (NN) 2 7i 
contact operator. Shown are d = 3 (red solid line), d = (black dashed line), and d = — 3 (blue dot-dashed 
line). The data are from Refs. [33-37]. 



Note, however, that the relatively large theoretical uncertainty of about 2M n /m^ « 30% for the 
cross section requires a carefull study of higher order effects. 

3. p-wave pion production 

Diagrams that contribute to p-wave pion production up to NNLO in the modified power count- 
ing are shown in Fig]3[ In particular, at NNLO there are subleading rescattering and direct pion 
production operators as well as the (AW) 2 71 contact term. Notice that it is the same contact term 
that also contributes to the three-nucleon force [8,23], to the processes yd — > nNN [24,25] and 
nd — ► yAW [26,27] as well as to weak reactions such as, e.g., tritium beta decay and proton-proton 
(pp) fusion [28,29]. Therefore it provides an important connection between different low-energy 
reactions. It is getting even more intriguing once one realizes that this operator appears in the above 
reactions in very different kinematics, ranging from very low energies for both incoming and out- 
going AW pairs in pd scattering and the weak reactions up to relatively high initial energies for the 
AW induced pion production. As a part of this connection in Ref. [30] it was shown that both the 
3 H and 3 He binding energies and the triton /3 -decay can be described with the same contact term. 
However, an apparent discrepancy between the strength of the contact term needed in pp — > pn% + 
and in pp — ► de + v e was reported in Ref. [31]. If the latter observation were true, it would question 
the applicability of chiral EFT to the reactions NN — ► NNn. To better understand the discrepancy 
reported in Ref. [31], in the recent paper [32] we simultaneously analyzed different pion produc- 
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Figure 5: Results for d 2 a / d£l K dM 2 p (left panel) and A y (right panel) for pn — > pp( l So)n ■ Shown are the 
results for ii = 3 (red solid line), d = (black dashed line) and d — —3 (blue dot-dashed line). The data are 
from TRIUMF [38, 39] (black squares) and from PSI [40] (blue circles) . 



tion channels. In particular, we calculated the /?-wave amplitudes for the reactions pn — > ppn~ , 
pp — > pnn + , and jop — ► c/7T + . Note that even in these channels the contact term occurs in entirely 
different dynamical regimes. For the first channel p-wave pion is produced along with the slowly 
moving protons in the l So final state whereas for the other two channels the l So pp state is to be 
evaluated at the relatively large initial momentum. In practice, we varied the value of the low- 
energy constant (LEC) d, which represents the strength of the contact operator, in such a way to get 
the best simultaneous qualitative description of all channels of AW — ► NNn. It should be stressed, 
however, that the value of d depends on the AW interaction employed and on the method used to 
regularize the overlap integrals. It therefore does not make much sense to compare values for d as 
found in different calculations. Instead one should compare results on the level of observables and 
this is what we do below (see also Ref. [32]). In Fig.|| we compare our results for various values 
of d with the experimental data for the analyzing power for the reaction pp — > dn + . We find that 
the data prefer a positive value of d of about 3. A similar pattern can be observed in the reaction 
pn — > ppK~ as illustrated in Fig. g. Again the data show a clear preference of the positive value 
for LEC d - our fit resulted in d = 3 for the best value. This channel, however, has been measured 
at TRIUMF at relatively large excess energy (17 = 0.66) where the conclusion may be spoiled due 
to the onset of pion d-waves. Fortunately, a new measurement for the same observables at lower 
energies is currently ongoing at COSY [41] which will soon allow a quantitative extraction of the 
value of the LEC d. We now turn to the reaction pp — ► pn7i + - this channel was used in the anal- 
ysis of Ref. [31]. The reaction pp — ► pnn + should have, in principle, the same information on 
the LEC d as contained in the deuteron channel. However, it is much more difficult to extract the 
pertinent information unambiguously from this reaction. In particular, the final NN-system might 
be not only in S- but also in P-wave both for isospin-zero and for isospin-one NN states. At the 
energies considered in the experimental investigation, v\ =0.22, 0.42, and 0.5, the Pp amplitudes 
may contribute significantly [43-45]. In the current study these states are disregarded. The results 
of our calculation for the magnitude A2 are given 2 in the left panel of Fig. 6 . One finds again 

2 The coefficients A/ are related to the unpolarized differential cross section via ^ = Aq + A2P2W , with Pi(x) 
being the second Legendre polynomial 
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Figure 6: Results for the mag- 
nitude A2 (left panel) and the 
partial wave amplitude ao(T i f]Ib) 
representing the relevant transition 
1 Sq — > 3 Sip (right panel) for pp — > 
pn7T + for different values of the 
LEC ii. The notation is the same 
as in Fig. 5, gray band corresponds 
to the results of Ref. [31]. The data 
are from Ref. [42]. 

that the positive LEC d ~ 3 seems to be preferred. Thus we conclude that all reaction channels of 
AW — > AW7T can be described simultaneously with the same value of the LEC d. Coming back to 
the problem raised in Ref. [31] it should be pointed out that the results of this work were not com- 
pared directly to the observables in pp — > pnn + . Instead, they were compared to the results of the 
partial wave analysis (PWA) performed in Ref. [42], as demonstrated in the right panel of Fig. 6. 
It is based on this discrepancy between data and theory it was concluded in Ref. [31] about the 
failure of simultaneous description of the weak processes and AW — > NNn. However, the partial 
wave analysis of Ref. [42] seems to be not correct. Here we refer the interested reader to Ref. [32] 
where the drawbacks of this PWA are discussed in detail. To illustrate the problems of the PWA 
in the right panel of Fig. 6 we also show the results of our calculation for the relevant partial wave 
ao which corresponds to the transition ^0 — > 3 S\p where the contact term acts. Clearly, although 
all data presented in Ref. [42] are in a good agreement with our calculation (see left panel in Fig. 6 
and also Ref. [32] for more details), the partial wave amplitude ao is not at all described. Thus, we 
think that the problem with the simultaneous description of pp — ► de + v e and pp — ► pnlt + , raised 
in Ref. [31], is due to the drawbacks of the partial wave analysis of Ref. [42]. 

4. CSB effects in pn -> dn° 

Recently, experimental evidence for CSB was found in reactions involving the production of 
neutral pions. At IUCF non-zero values for the dd — > an cross section were established [46]. 
At TRIUMF a forward-backward asymmetry of the differential cross section for pn — > dn° was 
reported which amounts to Aft, = [17.2 ± 8(stat.) ±5.5(sys.)] x 10~ 4 [12]. In a charge symmetric 
world the initial pn pair would consist of identical nucleons in a pure isospin one state. Thus the 
apparent forward-backward asymmetry is due to charge symmetry breaking. 

The neutron-proton mass difference is due to strong and electromagnetic interactions [47], 
i.e. SniN = m n — m p = Smp + 8m e ™. It was stressed and exploited in Ref. [11] that the strength 
of the rescattering CSB operator at LO in pn — > dn° (see Fig. 0(a)) is proportional to a different 
combination of Sniff and 8nf™ ( see also [48,49] for related works). Thus the analysis of CSB 
effects in pn — ► dn° should allow to determine these important quantities individually. It was, 
however, quite surprising to find that, using the values for 8mfj r and 5m™ from Ref. [47], the 
leading order calculation of the forward-backward asymmetry [11] over-predicted the experimental 
value by about a factor of 3 — a consistent description would call for an agreement with data within 
the theoretical uncertainty of 15% for this kind of calculation. The evaluation of certain higher 
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Figure 7: Leading order diagrams for the isospin violating s- 
wave amplitudes of pn — > dn°. Diagram (a) corresponds to 
isospin violation in the nN scattering vertex explicitly whereas 
diagram (b) indicates an isospin-violating contribution due to the 
neutron-proton mass difference in conjunction with the time- 
dependent Weinberg-Tomozawa operator. 

order corrections performed in Ref. [11] and in a recent study [50] did not change the situation 
sufficiently — the significant overestimation of the data persisted. In the recent work [51] we have 
shown that there is one more rescattering operator that contributes at LO (see diagram (b) in Fig. ^|). 
In full analogy to isospin conserving s-wave pion production, the idea was based on the fact that the 
energy-dependent WT nN vertex acquires an additional contribution proportional to 8m^ as soon 
as we distinguish between the proton and the neutron. We evaluated this new LO operator and 
recalculated the forward-backward assymetry atLO. It should be pointed out at this stage that Afb 
at LO is proportional to the interference of the s-wave pion CSB amplitude at LO and the p-wave 
pion isospin conserving amplitude. The latter is calculated up to NNLO, as discussed in Sec.|3], and 
exhibits very good description of data, which is a necessary pre-requisite for studying CSB effects. 
The complete LO calculation gives [51] 

<° = (11.5±3.5)xl0- 4 ^| (4.1) 

which agrees nicely with the experimental data if one uses the value of 8mf v from Ref. [47]. We 
may also use (Eq. [Ob to extract 8mff from the above expression using the data, which yields 

5m|f = (1.5±0.8 (exp.)±0.5 (th.)) MeV. (4.2) 

This result reveals a very good agreement with the one based on the Cottingham sum rule [47], 
Sniff = 2.0 ±0.3 MeV, and with a recent determination of the same quantity using lattice QCD [52], 
Smff = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV. 
5. Summary 

We surveyed the recent developments for NN — > NNn. We showed, in particular, (see Ref. 
[14, 16]) that the s-wave pion production amplitudes calculated up to NLO for pp — > dn + provide a 
good qualitative understanding of the pion dynamics. However, the relatively large theoretical un- 
certainty of about 2M n /mN « 30% for the cross section requires a computation of loops at NNLO. 
The latter are also absolutely necessary for pp — ► ppn , see Ref. [53] for the first results in this 
direction. Recently, we have studied /j-wave pion production up to NNLO [32]. In particular, we 
showed that it is possible to describe simultaneously the p-wave amplitudes in the pn — > pp%~ , 
PP —> pnn + , pp — > dn + channels by adjusting a single low-energy constant accompanying the 
short-range (NN) 2 n operator available at NNLO. We also demonstrated that the problem with the 
simultaneous description of the weak proton-proton fusion process and pp — ► pnn + , reported in 
Ref. [31], is most probably due to the drawbacks of the partial wave analysis of Ref. [42] used in 
Ref. [31]. Based on good understanding of the pion production mechanisms in the isospin con- 
serving case we studied charge symmetry breaking effects in pn — > dn°. We performed a complete 
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calculation of CSB effects at LO and extracted the strong contribution to the neutron-proton mass 
difference from this analysis. The value obtained, 8mff = (1.5 ±0.8 (exp.) ±0.5 (th.)) MeV, is 
consistent with the result based on the Cottingham sum rule and with the recent lattice calculations. 
At present the uncertainty in this results is dominated by the experimental error bars - an improve- 
ment on this side would be very important. On the other hand, a calculation of higher order effects 
is also called for to confirm the theoretical uncertainty estimate. 
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